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Abstract. The recovery of classical non-linear and chaotic dynamics from quantum 
systems has long been a subject of interest. Furthermore, recent work indicates that 
quantum chaos may well be significant in quantum information processing. In this 
paper we discuss the quantum to classical crossover of a superconducting quantum 
interference device (SQUID) ring. Such devices comprise of a thick superconducting 
loop enclosing a Josephson weak link and are currently strong candidates for many 
applications in quantum technologies. The weak link brings with it a non-linearity such 
that semi-classical models of this system can exhibit non-linear and chaotic dynamics. 
For many similar systems an application of the correspondence principle together with 
the inclusion of environmental degrees of freedom through a quantum trajectories 
approach can be used to effectively recover classical dynamics. Here we show (i) that 
the standard expression of the correspondence principle is incompatible with the ring 
Hamiltonian and we present a more pragmatic and general expression which finds 
application here and (ii) that practical limitations to circuit parameters of the SQUID 
ring prevent arbitrarily accurate recovery of classical non-linear dynamics. 
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Figure 1. Schematic of a simple SQUID comprising of a superconducting ring (of 
inductance L) enclosing a Josephson junction (with capacitance C). 

1. Introduction 

Quantum mechanics, like all physical theories, must satisfy the correspondence principle. 
That is, in the correct limit every physical theory must reproduce the predictions 
of classical mechanics. Perhaps because of its many counter-intuitive predictions, 
demonstrating that the correspondence principle holds for quantum mechanics has 
long been a subject of interest. One particular class of systems that has undergone 
extensive investigation comprise those that exhibit non-linear or chaotic dynamics. The 
interest here arising from the fact that the Schrodinger equation is strictly linear. This 
raises the question of how can non-linear phenomena emerge from such an equation 
for the evolution of the quantum state? A rather effective answer, that works for both 
dissipative [H [21 [Sj HI [5l [6] and Hamiltonian [7] systems is obtained when one introduces 
the effects of coupling to environmental degrees of freedom and simultaneously invokes 
the correspondence principle. Traditionally, studies of this kind have focused on the 
measurement problem and foundations of quantum mechanics. However, recent work, 
that builds on earlier work such as [HIE], has indicated that understanding such problems 
will be important in the development of new quantum technologies [TOl [HI [121 [13] . 
Indeed, a recent study has shown that a superconducting quantum interference device 
(SQUID) based on three Josephson junctions designed for the Josephson flux qubit may 
be used as a device for the study of quantum chaos [H]. Such work formed much of our 
motivation to study the recovery of non-linear trajectories for a related device that is a 
strong candidate for application in quantum information processing, namely the SQUID 
ring. In this work we consider a SQUID that takes the form of a thick superconducting 
ring enclosing a Josephson weak link as illustrated in figure [1] 

We note that while many textbooks present the correspondence principle as if 
it was a well understood and closed problem, the reality of the situation is far from 
this. Indeed, as we have already stated, in this paper we refine a large body of work 
based on trajectories methods that has seen a great deal of activity over the last 
fifteen years or so. However, not all systems, even when modelled classically, yield 
to a trajectory level analysis. For such systems a statistical treatment may be the 
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only possible approach. Hence, expressions of the correspondence principle fall into two 
categories namely the strong form (trajectory level) and the weak form (statistical and 
expressed in terms of distributions). While this paper is concerned with the strong 
form of the correspondence principle we take this opportunity to review an example 
of recent advances in understanding the quantum classical transition (QCT) in terms 
of the weak form presented by Greenbaum et. al. in [15]. The authors considered 
the quantum to classical transition of the chaotic Duffing oscillator. The classical 
dynamics is described though the evolution of a phase space distribution function that 
is the solution to the Liouville equation. This is compared with the dynamics of the 
Wigner function for the corresponding master equation. The authors show that the 
quantum to classical "transition occurs due to the dual role of noise in regularising 
the semiclassical Wigner function and averaging over fine structures in classical phase 
space" [IS]. Hence, correspondence was said to be achieved when cross-sectional slices 
of the Wigner function and the classical distribution function became, to all intents and 
purposes, indistinguishable. The authors considered the QCT in terms of scaling h - 
however we believe that their results are more general and that for non-linear systems 
the weak form of the QCT is achieved whenever the classical distribution function 
and the Wigner function become indistinguishable. Although their work is of great 
interest it does concentrate on the weak form of the correspondence principle. Finally, 
we observe that Greenbaum et. al. note in their work the conditions for weak QCT 
may need modification for systems with more than one dimension. Our experience 
modelling coupled quantum systems (especially modelhng atom-field collapse and revival 
phenomena and the QCT to Rabi oscillations [16j) indicates to us that expressions of 
the strong correspondence principle of the kind that we present in this work may well 
be valid for composite systems. However, this is a subject for extensive further study 
and we will not discuss it further here. 

We also take this opportunity to provide a brief overview of some other related 
work on the quantum to classical transition that we have not previously mentioned. D. 
Delande on p. 665 of [IT] and Friedrich and Wintgen in [18] discuss the emergence of 
chaos in the hydrogen atom. This is done within the context of scaling Planck's constant 
with respect to the classical action so that "at small (absolute) values of the energy the 
quantum spectrum becomes very dense, while it becomes sparse for larger field strengths 
and binding energies" [18j. Hence, the classical limit is said to be reached when the 
spectrum is dense and, hence, the effective value of Planck's constant is small. This is 
somewhat different from our approach where we scale the Hamiltonian (and therefore 
the spectrum) itself and also impose a localisation condition on the state vector. We 
also note that the work of M. Berry and B. V. Chirikov that are also in [T7j may well be 
of interest to the reader. There are a number of other important works, which mainly 
focus on solutions to the master equation and weak forms of the QCT, that nevertheless 
should be mentioned. The QCT for standard map and the quantum ratchet have been 
studied as dissipative driven systems in [19] and [20] respectively. We note that while [20] 
is of relevance to weak expressions of the QCT a trajectories approach was used as an 
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aid to computing solutions of the master equation. Finally, in [21] The Duffing oscillator 
is studied and shown to have a positive Lyapunov even when the system is far from its 
classical limit. Whilst this study focused on a master equation approach a trajectories 
technique was used in computing the Lyapunov exponents thus implying the emergence 
of chaos in terms of a strong form of the correspondence principle. 

SQUID rings, as with other superconducting devices, are unusual in so far as 
it is quite straightforward to perform experiments on them as a single macroscopic 
quantum object. Furthermore, they exhibit well known and understood semiclassical 
behaviour. Hence, as as individual trajectories are accessible, these devices form an 
ideal testing ground for better understanding the strong forms of the correspondence 
principle. Hence, this is the focus of this work. 

The correspondence principle in quantum mechanics is usually expressed in the 
form: 

"For those quantum systems with a classical analogue, as Planck's constant 
becomes vanishingly small the expectation values of observables behave like their 
classical counterparts" [223. 

However, and perhaps surprisingly, due to the ring's dependence on Planck's constant we 
find that this form of the correspondence principle cannot be applied here. Explicitly, we 
find that a problem arises when one considers the Josephson coupling term that is flux 
quantum ($0 = h/2e) periodic in external magnetic flux. The periodicity arises from a 
topological quantisation condition reached through an analysis of the Ginzburg-Landau 
Cooper pair wave function of the ring. This dependency of the flux quantum on h is 
highlighted by experiments that leverage this relationship in order to determine Planck's 
constant (see, for example, [23] )• Hence, we arrive at an interesting situation. In order 
to apply the correspondence principle we would normally wish to shrink h arbitrarily. 
However, we cannot do this without also changing $0? which is a topological, rather 
than canonical, quantisation condition. For later reference it is worth presenting a clear 
statement of the problem: 

There is a semiclassical description for the dynamics of the SQUID ring that 
is known as the resistively shunted junction (RSJ) model. It can predict, in 
some circumstances, results in excellent agreement with experiment. However, 
despite the Josephson current term containing an explicit dependence on 
Planck's constant, the model consists of equations of motion in terms of a 
classical co-ordinate (the flux). Hence, the solution to these equations of 
motion are classical trajectories. Our problem is then: From a quantum 
mechanical model of the ring, for a given solution to the RSJ model, how do we 
reproduce corresponding classical like trajectories in terms of the expectation 
values of observables? 

Instead of the usual expression of the correspondence principle we propose a more 
pragmatic and general statement which considers that Planck's constant remains fixed 
and that correspondence must be achieved through varying real circuit parameters. This 
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is expressed as a condition of large dynamics and localisation of the state vector when 
compared with a Planck cell in phase space. The localisation condition is key for systems 
whose solutions to the Schrodinger equation delocalise over time, such as the SQUID 
ring, as this requires the introduction of an additional localising effect. We note that 
unlike many other systems we find that practical limits on circuit parameters prevent 
us from obtaining an arbitrarily accurate correspondence between the quantum and 
classical dynamics. A feature that may well be of significance in quantum information 
processing applications of superconducting circuits. 

In this paper we start our discussion by considering the semiclassical resistively 
shunted junction (RSJ) model of the SQUID ring. This model is traditionally used, 
with confidence, when the ring capacitance is in the region 10"^^ to 10~^^F [21]. We 
then choose a typical set of circuit parameters and apply a sinusoidal external flux so 
that strongly non-linear dynamics are obtained. We show the Poincare section of the 
associated attractor. In doing this we follow the strategy of much of the previous work 
done on recovering classical like trajectories of non-linear systems [U, [31 HJ El E]. We 
then present the quantum description of this device. Obtaining classical like Poincare 
section's from the quantum model is achieved following our reformulated version of 
the correspondence principle. Here, localisation of the wavefunction arises through 
introducing environmental degrees of freedom via a quantum trajectories method, in 
this case, quantum state diffusion (QSD). This is an unravelling of the master equation 
that corresponds to a unit-efficiency heterodyne measurement (or ambi-quadrature 
homodyne detection) on the environmental degrees of freedom [25] (for an excellent 
introduction see [26], [27]). It is via such an open systems approach that we obtain both 
localisation of and non-linear dependence on the state vector from which we expect 
to be able to recover non-linear dynamics. Finally, we show that the quantum jumps 
unravelling of the master equation can be used to produce equivalent dynamics. This 
indicates that our results are not restricted to one environmental measurement process 
and may be quite general in nature. 

2. The RSJ model 

The SQUID ring is a device that exhibits macroscopic quantum degrees of freedom. 
However, under certain circumstances, the semiclassical resistively shunted junction 
(RSJ) model is often applied with a great deal of success. Although the RSJ model 
contains a term representing the purely quantum phenomena of tunnelling of the 
supercurrent through the Josephson junction it takes the form of a differential equation 
for a classical particle. That is, this model describes the evolution of classical trajectories 
and not state vectors. As with all quantum systems, the recovery of such classical 
trajectories should arise naturally from the application of the correspondence principle. 
The RSJ equation of motion for the magnetic flux, $, within a driven SQUID ring 
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where C, Jc, L and R are, respectively, the external flux bias, capacitance and critical 
current of the weak link, the inductance of the ring and the resistance representing the 
effects of environmental decoherence. The drive amplitude and frequency are Id and Ud 
respectively. Finally, $o = h/2e is the flux quantum. 

We make the following definitions: c^o = 1/V LC, r = c^o^; — $^.)/$o, 

(fi^ = $x/*o, /5 = 27rLIJ%, u = ujd/ooo, (fid = IdL/% and C = 1/2uqRC = 1/2Q where 



the quality factor Q = R^JC/L. We can then rewrite ([T]) in the standard, universal 
oscillator like, form 

^ + 2C^ + + sin [27r + (^,)] = sin {ujt) (2) 

We note now a property of this model that will become important later when we consider 
recovering classical like trajectories from quantum mechanical description of SQUID 
rings. That is, if we make the following transformation 

C ^ aC and L ^ hL (3) 

then we can keep the dynamics invariant by making the following changes at the same 
time: 

R \Jb/aR, Id Id/ and uJd — > oJd/ v^. (4) 

Even though non-linear and chaotic systems can exhibit extremely complex 
behaviour the Poincare section, comprising of points taken from the phase portrait once 
a drive period, characterises the dynamics. In figure [2] we show an example Poincare 
section for a SQUID ring with circuit parameters C = 1 x 10~^^F, L = 3 x 10~^°H, 
R = lOOQ, P = 2, Ud = ujq and Id = 0.9 /lA. We note that we have biased the ring at 
the half flux quantum, = 0.5$0; so that the potential approximates a double well. 

3. The quantum model and correspondence 

The quantum description of the SQUID ring is obtained when one considers the following 
standard Hamiltonian [28] : 

2^ 2L 2e V "^0 / 

which is consistent with the equation of motion ([1]). Substituting $o = h/2e in order to 
explicitly show the Hamiltonians dependence on h yields, 

H^9lJlZ^.H,„J?^) (6) 

2C 2L 2e \ h J ^ ^ 

where the magnetic flux threading the ring, $, and the total charge across the weak link 
Q take on the roles of conjugate variables for the system with the imposed commutation 
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if (magnetic flux in units of ^o) 



Figure 2. Poincare section obtained from the RSJ dynamics for a SQUID ring with 
circuit parameters C = 1 x IQ-^^F, i = 3 x lO-^^H, R = lOOQ, /3 = 2, = wq, 
Id = 0.9 /iA and = 0.5$o- 



relation Q = ih. For clarity we will use the circumflex to denote operator status 
throughout this paper. Here is the external applied magnetic flux and incorporates 
the drive term for the ring. 

We note that from the Hamiltonian one might naively consider that in the limit 
^ — > that the Josephson junction coefficient ?ilc/2e becomes neghgible and the system 
reduces to a harmonic oscillator, where both classical and quantum models could them 
be made to agree. However, as previously noted, the RSJ model produces a very good 
match to experiment for situations where hlc/2e (and therefore h) is non-negligible. It 
is a correspondence with such predictions we are interested in recovering. 

It is usual, when working with SQUID rings, to normalise flux and charge operators 
in terms of $o and 2e respectively. However, as $o is dependent on Planck's constant 
and 2e is not we will not adopt this convention here as doing so obscures the role of h. 
Instead, we define dimensionless fiux and charge operators in the manner usual for the 
simple harmonic oscillator: x — ^JCuJoJh^ and p — ^Jl/hCuJaQ. These quantities 
give us a clear framework from which to discuss the correspondence principle. To 
see why this is the case, let us consider the fact that the commutator |^$, Q] = ih, 
which implies the uncertainty relation A$AQ > h/2. This sets a minimum area of 
uncertainty in phase space, the Planck cell, of h/2. Hence, in the conventional view 
of the correspondence principle, letting ^ — is equivalent to making the Planck cell 
vanishingly small. This implies that as /i — > states localised to a Planck cell will look 
more and more like point particles and effects such as quantum noise fluctuations and 
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tunnelling will become less and less significant. In dimensionless units the uncertainty 
relation becomes AxAp > 1/2. We will, for future reference, associate a characteristic 
length with each quadrature of Sx = 5p = 1/V2, i.e. the length of one side of a square 
of area equal to a Planck cell. Hence in this representation the size of the Planck 
cell is fixed. It is now the system that scales, in this case relative, to 6x and 6p. For 
example consider the Duffing oscillator, a non linear system with a double well potential. 
When described in dimensionless units the effect of reducing h is to deepen and broaden 
the wells in the potential. Hence, increasing the scale of the systems dynamics with 
respect to the now fixed characteristic length 6x. It is this change of scale together 
with the introduction of an environment that localises the system to a Planck cell that 
has enabled previous applications of the correspondence principle in terms of h to be so 
successful m El [am [HIE]. 

Returning to our discussion of the ring Hamiltonian we note that the first two 
terms are the (shifted) simple harmonic oscillator potential. Hence, normalising H 
by hujQ would imply that the curvature, and therefore spectrum, of such an harmonic 
oscillator would remain unchanged if h is scaled. Thus, in this system of dimensionless 
units and normalised energy we can use Sx, 6p and the background parabolic curvature 
as a constant reference by which to gauge the effect of changing h. This is apparent 
when we explicitly give the dimensionless Hamiltonian: H' = H /Tiuq we find that 



H = h cos 

2 2 2euon 




1/2 
X 



(7) 



Where x{t) represents the normalised external bias fiux together with a drive term 
corresponding to the one used in the RSJ model. We define, for our later convenience, 

the multiplying factor of x within the cosine to be i7 = [Ae^ /h)J{L/C) 

It is now that we can see why the SQUID ring does not scale appropriately with h. 
Our discussion of this point makes use of an on-line animation that accompanies this 
paper. The reader is strongly advised to make use of it as it will greatly clarify the 
behaviour of the ring as Planck's constant is reduced. 

In the animation we show the potential energy of the ring together with the energy 
eigenvalues. The eigenvalues are plotted only where the probability density function of 
the eigenstates are greater than 0.05. Thus we can identify some salient features of the 
eigenfunctions, such as nodes, without having so much information available that the 
animation is unclear. Planck's constant is reduced by one percent per frame until it is 
scaled to approximately 2 x 10^^ of its original value. We note that this movie uses the 
dimensionless units described above. Hence the characteristic length 5x = and 
the background parabolic curvature associated with the harmonic oscillator part of the 
Hamiltonian can be used as a reference in fiux and energy respectively. 

At the beginning of the animation we see a double well potential. These wells 
initially deepen and move closer together with more energy levels being available in 
each well. About twelve seconds into the animation the wells approach their maximum 
depth but nevertheless continue to narrow. The result of this is that the number of 
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energy levels in each well begin to reduce. Shortly after this point, as Vt increases, 
more wells close in from the outside of the potential. As the animation continues, more 
and more wells appear in this plot and there are ever fewer energy levels contained 
in each well (with the zero point energies of each well rising significantly). We also 
note that there is some tunnelling near the top of the wells. The trend continues 
until approximately one minute into the animation when the lowest energy eigenstates 
start to delocalise over a number of wells beyond the central double well. This is due 
to significant tunnelling under the barrier as the distances between the wells becomes 
smaller. This trend continues and at the end of the animation we see that the spectrum 
seems to comprise of two parts - that within the wells and an almost harmonic oscillator 
type component due to the parabolic background. 

Clearly the potential is scaling the wrong way as — > and we are moving away 
from a situation where we can recover the type of dynamics seen in Fig. [21 Indeed, by 
looking at the way the spectrum is changing it appears that the system is becoming 
manifestly more, rather than less, quantum mechanical. Furthermore, if we were to now 
introduce a localising environment following [H [21 El IH El E] then the noise fluctuations, 
being significant compared with a Planck cell, would simply wash out any details. In this 
way recovery of a Poincare section as shown in figure [2] will become impossible. Hence, 
it is evident that the usual expression of the correspondence principle is incompatible 
with our objective of recovering RSJ type trajectories for the SQUID ring. 

As the notion of correspondence is of fundamental importance we now try to recover 
the situation and propose a slightly different statement of the correspondence principle 
from the one already given: 

Consider h fixed (it is) and scale the Hamiltonian, in such a way as to 
preserve the form of the classical phase space, so that when compared with the 
minimum area ^/2 in phase space: 

(a) the relative motion of the expectation values of any observable 
(generalised) co-ordinate (and hence the associated classical action) 
become large and 

(b) the state vector is localised (in any representation)|j]. 

Then, under these circumstances, expectation values of operators will behave 
like their classical counterparts. 

X We note that as this item is phrased in terms of an area in phase space that the locahsation condition 
must apply to aU quadratures. Hence, in order to aid conceptuaUsation of this condition it may be 
helpful to consider it in terms of the Wigner function. Here, we would require this distribution to 
approximately cover, and be localised to, the minimum area h/2. Non-negativity of the Wigner function 
is also sometimes stipulated as a condition for the QCT. However, fluctuation terms from interaction 
with the environment may generate short lived negative sections of the Wigner function. Here we will 
not impose non-negativity of the Wigner function as a condition for the QCT as these interference 
terms will, in the correspondence limit, be local and small scale. Therefore, they will have a negligible 
effect on the systems large scale dynamics. 
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We note that for many Hamiltonians that part (a) here is mathematically equivalent to 
the original statement of the correspondence principle. That is, as already discussed with 
reference to the Duffing oscillator, reducing h is the same as scaling the Hamiltonian. 
However, it is not the only way to scale the Hamiltonian as we can also change the 
systems parameters. 

We have seen that the statement (a) is not the same as reducing h for the SQUID 
ring. In order to satisfy our new statement of the correspondence principle we see that 
we must instead scale the Hamiltonian by adjusting circuit parameters. This must be 
done in such a way so as to preserve the form of the dynamics of the classical phase 
portrait (i.e. the phase portrait may be scaled relative to a Planck cell but not distorted 
in any other way). It is now that we see the value of writing the RSJ equation of motion 
in the dimensionless form ([2]). That is, if we change the capacitance and/or inductance 
using ([3]) and simultaneously make the changes (jll) that the equation of motion ([2]), and 
the corresponding phase portrait, remains unchanged. However, we can also see that in 
terms of the dimensionless flux, (x), and charge, (p), that the Hamiltonian [7] will scale 
appropriately for the QCT for increasingly large values of capacitance. Hence, we have 
found a mechanism, via the parameter a, for scaling the ring Hamiltonian with respect 
to the plank cell that avoids the problem we encountered when letting h tend to zero. 

We now have an expression of the correspondence principle which encompasses a 
larger class of systems than under the original expression. Unfortunately, this alone 
is insufficient to recover classical trajectories. For systems such as the SQUID ring 
solutions to the Schrodinger equation will delocalise to cover much of the possible phase 
space and expectation values will not behave like their classical counterparts. This can 
be seen clearly if one looks at phase space representations of the wave function [29l [T5] . 
In order to overcome this difficulty we have introduced the extra constraint (b) of 
localisation. In this work we localise the wavefunction by treating the SQUID ring as an 
open quantum system where localisation arises due to the introduction of environmental 
degrees of freedom. 

In an ideal world we would have total control of the parameters of the SQUID ring. 
Hence, application of the above expression of the correspondence principle could be 
implemented to arbitrary precision. However, the ring is a real physical device and the 
above expression of the correspondence principle will be limited by achievable circuit 
parameters. Typically C takes values between 10~^^F and 10~^^F with the RSJ model 
being applied for C = lO^^^F and larger. These constraints together with the fact that 
ri scales as {L/Cy^'^ very much restricts the scalability of this Hamiltonian. That is, 
there is a practical limit to how classical this device can become. 

In any statement of the correspondence principle we expect the expectation values 
of the observables to behave the same as those of the classical counterpart. By using 
the Hamiltonian of Eq. ([71) and introducing dissipation as done in the next section we 
can not produce results that are in agreement with those from the classical equation 
of motion ([2]). This is due to the fact that this method of introducing Ohmic damping 
does not bring with it the frequency shift that arises through the damping term in the 
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classical dynamics. We can easily resolve this problem addition of an extra term to the 
Hamiltonian [22l [26] . That is, the Hamiltonian 

2 + 2^ ~ 2^ ^ ^ + \{P^ + ^P) (8) 

is one we can use to generate trajectories that will be comparable to those predicted by 
the RSJ model. 



4. QSD and the recovery of classical dynamics 

We now introduce the QSD description for the evolution of the state vector that is 
given by the Ito increment equation [301 

= - ^H' 1^) dt 



+ E 



+ E [4 



0) dt 

di (9) 



where the first term on the right hand side of this equation provides Schrodinger 
evolution and the second and third terms describe the decohering effects of the 
environment on the evolution of the systems state vector. In this work we consider 
only one Lindblad operator L = y/2(d, where a is the oscillator lowering operator. The 
time increment is dt and the dC, are complex Wiener increments that satisfy rf^^ _ _ q 
and d^d^* = dt [30l[3l] (the over-bar denotes the average over stochastic processes). 

In figure [3] we show the computed Poincare section obtained for (x) and (p) 
using this quantum state diffusion model. Here the SQUID ring has the same circuit 
parameters as given in figure [H We can see that although we have not reproduced 
exactly the Poincare section of figure [2]the same structure appears to be indicated here. 

From our discussion above we would therefore expect to improve our resolution of 
this structure for values the scaling parameter a > 1 where a is defined in equation ([3]). 
We choose a = 100 so that we obtain a capacitance, C = 10~^^F, commensurate with 
the the upper practical limit of typical Josephson junctions. In figure H] we show the 
Poincare section for this choice of a. It is here that the fact that scales as C~^^^ 
becomes important. As we can clearly see we have not improved the resolution of detail 
to the extent that we would have liked. Previous studies of the Duffing oscillator find 
good correspondence for trajectories [32l [33] that cover approximately three times the 
area of phase space than exhibited in these dynamics. In order to attain this limit we 
would need a ring capacitance of the order of 10~^F. Even then some quantum noise 
would be evident and true convergence to the RSJ model will be achieved only for still 
greater, and more unrealistic, values of capacitance. 
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Figure 3. Poincare section obtained from the quantum state diffusion model for 
the SQUID ring with circuit parameters as given in figure [2l As (x) and (p) are 
dimensionless, we have indicated on top of this figure the magnetic flux in so that 
direct comparison can be made with figure [2l 



5. Quantum Jumps 

For the preceding results, we have used the QSD unravelhng of the master equation. 
This unravelhng makes assumptions about the underlying measurement process applied 
by the environment. For QSD this is unit efficiency heterodyne detection (or ambi- 
quadrature homodyne detection) [25]. Hence, one can only draw limited conclusions 
from a demonstration of the quantum to classical transition using a specific unravelling. 
In order to see if our observations of the correspondence limit for the SQUID ring may be 
general in nature we choose to reproduce the results of figure H] using another unravelling 
of the master equation - quantum jumps [Ml 135] . This model is very different from QSD 
as it is based on a discontinuous photon counting measurement process, rather than a 
diffusive continuous evolution. In quantum jumps the pure state stochastic increment 
equation is given by 

\dilj) = - ^H' lip) dt 



dt 
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Figure 4. Poincare section obtained from tlie quantum state diffusion model for 
the SQUID ring with circuit parameters changed from Fig. O by setting a = 100 
as introduced in equation ([3]). This corresponds to increasing the capacitance to 
C ~ lO^^^F. For clarity we note that this Poincare sections width is in terms of 
flux quantum approximately the same as that in the section presented in figures 1 and 
2. The change in scale, in terms of the dimensionless flux {{x)) and charge ((p)), for 
this plot has arisen through the scaling parameter a. 



+ E 



(10) 



where dNj is a Poissonian noise process satisfying dNjdNk = 6jkdNj, dNjdt = and 



dN, 



dt. That is jumps occur randomly at a rate that is determined by [LjLj 



Again the Lindblad operator is L = ^/2^d. 

If we were to hmit our discussion of the QCT of SQUID rings to QSD our 
results would be meaningful only when the environment operates as a continuous unit 
efficiency heterodyne detection (or ambi-quadrature homodyne detection); however, by 
applying another unravelling such as quantum jumps we can investigate if our results 
are measurement process dependent. We find that the Poincare sections shown in 
figures m and [5] generated using QSD and quantum jumps (respectively) are in very good 
agreement. While it is not possible to model the infinity of possible unravellings of the 
master equation; the fact that two very different models of environmental measurement 
processes produce such similar results is, to us, strongly indicative that our results are 
quite general in nature. 
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Figure 5. The Poincare section using identical parameters as those used to compute 
figure [4] where the evolution is now computed using the quantum jumps unravelling of 
the master equation instead of quantum state diffusion. 



6. Conclusion 

In this work we have presented a new, pragmatic, reformulation of the strong form of 
the correspondence principle. The key features being scaling the Hamiltonian (but 
not through changing h) and a localisation of the state vector. Explicitly, in this 
formulation it is required that the relative motion of the expectation values of observable 
co-ordinates become large and the state vector is localised when compared with the 
minimum area ^/2 in phase space. Then, under these circumstances, expectation values 
of operators will behave like their classical counterparts. We have used this expression 
of the correspondence principle to understand the QCT for SQUID rings. Here the 
classical limit is understood to be correspondence of the dynamics of expectation values 
of the quantum system with the semi classical RSJ model. To the best of our knowledge 
this is the first time that the QCT for SQUID's has been demonstrated within a general 
expression of the correspondence principle. Because we see a QCT for SQUID's it is our 
opinion that this work reinforces the point of view that the weak or strong formulations 
of the correspondence principle should not be phrased in terms of scaling h. 

When our statement of the correspondence principle is applied to the SQUID 
ring we achieve some success in recovering recognisable classical non-linear trajectories. 
However, this is somewhat mitigated by practical constraints on circuit parameters. 
To be precise the Josephson term in the Hamiltonian scales as the fourth root of the 
ratio of the inductance to the capacitance. It is this dependence that, in practice, 
prevents arbitrarily good correspondence between the predications of the RSJ model 
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and both the QSD and quantum jumps unravelhng of the master equation. Hence we 
note that quantum effects are therefore non-neghgiblc in the region of parameter space 
where semi-classical models would normally be used with confidence. One immediate 
conclusion that can be drawn from this observation is that the the RSJ model will be 
more effective in situations where any drive amplitude on the ring is large. A good 
example of such a situation are voltage dynamics experiments where, historically, the 
RSJ model has been of great utihty. Furthermore, this expression of the correspondence 
principle may well turn out to be of importance when considering Josephson junction 
devices for application within the field of quantum information processing. 
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